Quantum Hall fractions in rotating Bose-Einstein condensates 
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We study the Quantum Hall phases that appear in the dilute limit of rotating Bose-Einstein 
condensates. By exact diagonalization in a spherical geometry we obtain the ground-state and low- 
lying excited states of a small number of bosons as a function of the filling fraction v, ratio of the 
number of bosons to the number of vortices. We show the occurrence of the Jain principal sequence 
of incompressible liquids for v = 1/2, 2/3, 3/4, 4/3, 5/4 as well as the Pfaffian state for v = 1. The 
collective excitations are well described by a composite-fermion scheme. 
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Bose-Einstein condensates in dilute atomic gases offer 
a unique opportunity to investigate the physics of vortex 
matter when they undergo rotation^, ^]. Indeed, recent 
experiments j^, have observed the appearance of large 
vortex arrays at sufficient high angular velocity u>. In 
addition to this phase akin to the Abrikosov lattice of 
type-II superconductors, there is the possibility that at 
larger ui the lattice melts [^j and is replaced by a quan- 
tum Hall liquid. Consider a trap with strong confinement 
in the z direction such that the system is effectively two- 
dimensional (2D). Then if the rotation frequency is tuned 
to the characteristic frequency of the harmonic confining 
potential in the xy plane, the bosons feel only the Coriolis 
force and the system is equivalent to 2D charged bosons 
in a magnetic field, i.e. the conditions of the quantum 
Hall effect. In this regime, it has been pointed out H, |tJ 
that the celebrated Laughlin wavefunction is the exact 
ground state for the filling fraction v = 1/2 where v is 
the ratio of the number of bosons to the number of vor- 
tices. Some of the excitations above this ground state are 
quasiparticles with fractional statistics which may even- 
tually be probed by laser manipulations B). Investiga- 
tions by exact diagonalization have given evidence for 
even more exotic states of matter ]To|, pd| , some involving 
parafermionic wavefunctions introduced in the context of 
the fractional quantum Hall effect for fermions Jl2| . 

In this Letter we investigate the quantum Hall states 
of bosons as a function of the filling v by use of exact 
diagonalizations in the spherical geometry jyl |l4| . This 
allows to separate bulk from edge excitations. We show 
the appearance of the Bose analog of the Jain princi- 
pal sequence of fractions, v — ;^pj-,^zy- The excited 
states show collective modes well described by a compos- 
ite fermion picture in which there is binding of one vor- 
tex per boson. We obtain evidence for the Pfaffian state 
|jlO| , |TT[| at v = 1 by displaying its peculiar half-vortex 
excitations. For higher fillings, v > 3/2, we observe 
some states with properties of the Read-Rezayi (RR) 
parafermionic states. However they show no clear ten- 
dency to convergence to the thermodynamic limit. 

In the rotating frame Jig] , the Hamiltonian describing 



N bosons of mass m is given by : 
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where the xy trap frequency is ujq, the axial frequency is 
u) z and the angular velocity vector is ujz. In the ultra- 
cold atomic gases the interaction takes place through 
s-wave scattering only and is thus given by V(r) — 
/m)6^ 3 \r) where a s is the s-wave scattering 
length. For uj close to ljq, the physics is that of charge-e 
bosons in a magnetic field B = (2mu/e)z, correspond- 
ing to a magnetic length £ = y/h/ (2mw). There is then 
a 2D regime in which the boson wavefunction along the 
z-axis is the ground state of the harmonic oscillator and 
the interaction is now given by V 2D (r) — g£ 2 8^(r) (the 
vector r is 2D) 



with 



// 32irTiuja s /£ z , l z — ^Jfi/mu) z 
is the confinement length along z. The energy scale of 
the quantum Hall problem is thus set by g. 

We are thus led to study the quantum Hall effect of 
bosons interacting through a delta potential in the low- 
est Landau level (LL) pq . To study the vortex liquids 
that appear as a function of the filling factor v, we use the 
spherical geometry [fl4| , [l7) in which the bosons move on 
a sphere of radius R in the magnetic field of a monopolc 
B = hS/eR 2 at the center of the sphere, giving rise to 
2S+1 cyclotron orbits in the lowest LL. In the thermody- 
namic limit, the filling factor v is given by N/2S. How- 
ever for the incompressible liquids there is in general a 
finite shift in the relation between the number of particles 
and the flux, i.e. one has generally 25 = (l/v)N — X. 
If we have a guess for the ground state then one can 
evaluate the shift and check for its validity against nu- 
merical results. For example the bosonic Laughlin state 
for v = 1/2 on the sphere is realized for 2S — 2N — 2 by 
the wavefunction : 
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where the spinor coordinates (it, v) are given by : 

( Ui , Vi ) = (cos^/2e^/ 2 ,sin^/2e-^/ 2 ). (3) 

This is an exact zero-energy eigenstate of the present 
problem 0. We have conducted Lanczos diagonaliza- 
tions for various N and flux 2S to elucidate the nature of 
the incompressible liquid states. States can be labelled 
by their total angular momentum L, contrary to the pla- 
nar geometry where only the z component is conserved. 

Jain sequence. The signature of incompressible states 
is the presence of a L=0 singlet ground state separated 
by a clear gap from excited states. A typical spectrum 
is given in fig. la for v = 1/2. In the excited states we 
observe a well-defined collective mode which is gapped 
for all values of L. Clear signs of incompressible liq- 
uids are seen for the fractions v = 1/2, 2/3, 3/4, 4/3, 5/4 : 
some sizes are displayed in figs. lb-d. These fractions are 
the bosonic analog of the Jain sequence jl8| of fractional 
quantum Hall states. They are explained by a composite 
particle picture in which the composite fermion (CF) is 
a boson bound to one vortex. Then the integer quan- 
tum Hall effect with n filled CF LLs leads to a fraction 
at v = — tt. This state is realized on the sphere for 
2S = (n+ l/n)N — (n + 1) which is exactly what we ob- 
serve. The collective mode is then an exciton-like mode 
obtained by promoting one CF from the highest occu- 
pied LL to the next LL. On the sphere the maximal L is 
then given by L max = N / n + n — 1 in complete agree- 
ment with our results. We estimate the gap to these 
neutral excitations by finite-size scaling [|l9| : ~ 0.09g 
(v = 1/2), 0.05g (2/31 0.04g (3/4). The series of values 
for (N, 2S) is aliased with the sequence of fractions at 
v = Indeed, if (N, 2S) matches the fraction n/n+1 
then it also matches fraction p/p — 1 for p=N/n. Hence 
the same data set points to the presence of the fractions 
v = 4/3 and 5/4 (we do not have enough points to pro- 
vide gap estimates). The two-particle correlation func- 
tion g(r) is displayed in fig. 2 for some states in the CF 
sequence. For v = 1/2 it is essentially free of finite-size 
effects and shows the strong correlation hole characteris- 
tic of a Laughlin ground state. The state with v = 2/3 
no longer vanishes at the origin and hence has a nonzero 
ground state energy. The CF sequence also include the 
fractions v — 3/2 and v = 2. For these values we find 
families of incompressible states in the (N, 2S) plane but 
they show no sign of convergence toward the thermody- 
namic limit, neither in the ground state energies nor in 
the gap values. For these fractions, we have candidates 
possibly originating from the Read-Rezayi parafermionic 
wavefunctions (see below). 

For fillings less than 1/2, the spectrum has many zero- 
energy eigenstates that are the quasiholes of the Laughlin 
state eq.(^]). This is a special property of the delta func- 
tion interaction, the " quasielectron" being gapped. This 
obscures the appearance of fractions less than 1/2. If we 
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Figure 1: Energy spectrum for (a) 8 bosons at 2S=14 (u — 
1/2), (b) 8 bosons with 2S=9 (2/3), (c) 12 bosons at 2S=12 
(3/4), (d) 8 bosons at 2S=5 (4/5). Energies are in units of g 
and the horizontal axis is total angular momentum. 
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Figure 2: Two-particle correlation function g(r) as a function 
of great circle distance on the Haldane sphere in units of the 
magnetic length. The v = 1/2 curve is plotted for sizes N=8,9 
and for 2/3 N=9,10. 



change the interaction from pure delta by adding a pseu- 
dopotential V% H] in the next allowed partial wave for 
bosons, i.e. I = 2, then we find other states from the 
hierarchy, the strongest being v = 2/5. 
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Figure 3: Gaps for the Bose Jain sequence 1/2 and 2/3 as 
well as the Pfaffian state. 



Pfaffian state. The filling v = 1 corresponds to the 
absence of magnetic field acting upon the composite 
fermions. Previous studies p(J are indicative that pair- 
ing of the CFs takes place instead of a Fermi surface. 
An appealing wavefunction describing this phenomenon 
is the so-called Pfaffian state. On the sphere it can be 



Figure 4: (a) Spectrum at the Pfaffian matching condition 
N=12, 2S=10, with a collective mode above an isolated sin- 
glet ground state (b) with one extra flux quantum, two quasi- 
particles give rise to a degenerate set of states L=0,2,4,6. 



written as : 

#„=i = Pf{ — - — }TT( M ^ - u i v i)> ( 4 ) 



UiVj — UjV l 



i<3 



where Pf stands for the Pfaffian which is the antisym- 
metrized product of pair wavefunctions pl[ (a fermionic 
version of this state is a good candidate to describe 
the enigmatic v = 5/2 quantum Hall state). Calcula- 
tions of overlaps between the model Pfaffian wavefunc- 
tion Eq.(||)and the exact ground state suggest that it 
describes the physics of bosons at v = 1 in toroidal and 
disk geometry || (?], |). The bosonic Pfaffian state is 
realized on the sphere for 2S=N-2 for all N even. Our 
calculations lead to incompressible states at these spe- 
cial values for N=4,6,8,10,12. There is a clear gap that 
extrapolates smoothly to w 0.05g. It appears on fig. 4a 
for N=12. This state has charged excitations that are 
different from those of a Laughlin fluid. If we add or 
remove one flux quantum, then two quasiparticles are 
created, leading to a set of low-lying states with an alter- 
nate even-odd character : see fig. 4b. This is consistent 
with the spectrum for two identical particles with repul- 
sive interactions. This is observed for all accessible sizes. 
We consider this as a proof that the physics is different 
from that of the CF sequence and is the hallmark of the 
Pfaffian state The correlation function is shown in 
fig. 5. It has now a hump at the origin possibly due to the 
pairing of the CFs. The energy of this state is lowered 
by a correlation dip (instead of a hole) which appears at 
some characteristic radius of order one magnetic length. 
For larger separation, g(r) approaches 1 but with a char- 
acteristic length scale which is definitely larger than that 
occurring in fig. 2. The correlation length of Pfaffian is 
larger than for the Jain-like fractions. 

Read-Rezayi states. For larger fillings it has been sug- 
gested by Cooper et al. Q that fractions occur at v = k/2 
and are well described by the Read-Rezayi parafermionic 
wavefunctions. These functions involve clustering of k 
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Figure 5: The two-particle correlation function for the exact 
ground state at the Pfaffian point 2S=N-2. It has an extra 
concentration of bosons at the origin. 



particles and are a generalization of the Pfaffian which 
corresponds to simple pairing, i.e. k = 2. On the sphere 
they are r ealizedfor 2S = fiV-2. There are possible can- 
didates at v = 3/2 for N=6,9,12,15, v = 2 for N=8,12,16 
and the fraction v = 5/2 may be realized for N=15 and 
20 (but the gap is very small). Contrary to the samples 
belonging to the CF sequence (they have different fluxes 
since the shifts are different between the CF sequence and 
the RR states) we find no sign of convergence towards a 
thermodynamic limit. The gap is non-monotonous as a 
function of the size for v = 3/2, and for v = 2 finite-size 
effects are very large, preventing any sensible extrapola- 
tion. One possibility is that these states have very large 
correlation lengths and are not accommodated on our 
largest spheres. This is consistent with the fact that the 
correlation function shows very strong oscillations and no 
hint of incompressibility. 

Evidence for clustering of more particles comes from 
the correlation function where we see the same phe- 
nomenon as in the Pfaffian case. The hump at the origin 
is even more pronounced and surrounded by a correlation 
hole. The hump also increases with the filling albeit we 
cannot make quantitative statements. 

If we increase the number of bosons at fixed flux, 
then the spectrum becomes rotor-like : the levels lie on 
parabolas described by effective Hamiltonian 157 L 2 and 
correlation effects disappear. 

We have shown the appearance of the Jain princi- 
pal sequence of quantum hall fractions in rotating Bose- 
Einstein condensates. The composite fermion picture 
gives a successful account of the observed fractions as 
well as their collective mode excitations. The Pfaffian 
state is realized at v = 1 as seen from the special match- 



ing of flux and number of bosons as well as its half-flux 
quasiparticle. The gaps we estimate from our diagonal- 
izations are all of the order of Tiuja a /l z . 

We thank Yvan Castin and Jean Dalibard for numer- 
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